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(g) Pro cede et dispositif permettant a un modem de se synchroniser sur un transmetteur de donnees 
numeriques par voie hertzienne en presence de brouilleurs. 
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(57) Le procede consiste a estimer (9) la matrice 
de correlation ftxx des signaux recus du trans- 
metteur sur un ensemble de n capteurs de 
reception, a calculer (10) la matrice de correla- 
tion inverse & calculer (11) des vecteurs 
d'intercorrelation t Mi entre les signaux recus 
sur I'ensemble des n capteurs et un signal de 
replique d(k) connu, a calculer un critere de 
synchronisation (15) en effectuant les produits 
des vecteurs dMntercorrelation transposes 
conjugues ftii et de la matrice des vecteurs 
d'intercorrelation fxdi. et>a comparer (16) la 
valeur du critere obtenu a une valeur de seuil i\ 
determinee pour placer la synchronisation sur 
rechantillon pour lequel la valeur des criteres 
depasse la valeur de seuil j\. Pour estimer le 
decalage Doppler, la sequence de synchronisa- 
tion est' decoupee en tranches d'un nombre 
determine de symboles pour effectuer des cat- 
culs de transformation de Fourier rapide sur 
des vecteurs d'intercorrelation determines sur 
chaque tranche. 
Application : Transmission radio de donnees. 
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La presents invention cencerne un procede et un dispositif permettant 'aun modem de se synchroniser 
sur un transmetteur de donnees numeriques par voie hertzienne notammeht en presence de brouilleurs. 

Dans les systemes de transmissions de donnees numeriques par voie hertzienne la synchronisation est 
effectuee a partir de sequences de symboles connus du recepteur et intercalees entre les symboles d'infor- 
5 mation par correlation entre le signal recu et un signal de reference construit a partir des symboles connus. Le 
resultat est compare a un seuil pour decider de la synchronisation. Des systemes permettant la synchronisation 
d un modem a partir d'un seul capteur et d'un seul recepteur sont par exemple connus de I'article D.V. Sarwate et 
M.B. Pursely intitule "Cross correlation properties of pseudo random and related sequences" paru dans IEEE Vol. 
66, pp. 593-619, May 1980. Mais les techniques mises en oeuvre qui permettent en outre une estimation des ca- 
10 racteristiques du canal de transmission deviennent rapidement inef f icaces en presence de brouillage. 

Lorsque le rapport signal a bruit n'est pas suffisant pour la prise de synchronisation par une technique 
mono-capteur, des techniques visant a accroitrele rapport signal a bruit et a permettre la synchronisation sont 

utilisees. , ,. .. . . . 

Ces techniques peuvent etre regroupees en deux families. Une premiere famille vise a pretraiter le signal 
is recu avant d'effectuer la synchronisation du signal pretraite. La seconde famille vise a ameliorer le rapport 
signal a bruit et a effectuer la synchronisation conjointement. 
T Dans les techniques relevant de la premiere famille il faut egalement distinguer deux categories : 
' ,; '• ' '•: une premiere categorie met en oeuvre un pretraitement mono-capteur de type excision de brouilleurs 
;v ; ' : ' avaht d'effectuer la synchronTsatlbn- C*-tiaTteiir(erit fesfcbnnu-par exemple de I'artide de LB. MILSTEIN 
20 "Interference rejection Techniques in spread spectrum communications", paru dans Proc. IEEE, Vol. 76, 

pp 657-661 June 1988. Suivantcet article on cherche a predire la contribution du brouilleur, dont la 
bande est supposee plus petite que la bande du signal, sans predire celle du signal utile. Cela permet 
d'ameliorer les performances de la synchronisation en presence d'interferences bande etroite (type 
FSK) Mais le traitement reste ineff icace en presence d'un brouillage occupant toute la bande utilisee 
25 par le signal utile. Cela est en partie due au fait que les techniques de pretraitement mono-capteur ne 

permettent un traitement des signaux que dans un espace limite a I'axe des temps et a I'axe des fre- 
quences. \ 
- La seconde categorie concerne les techniques mettant en oeuvre un pr6traitement multi-capteurs de 
rejection, precedent la prise de synchronisation. L'idee consiste a mettre en oeuvre un f iltrage spatial, 
30 ' " discriminant le signal et les brouilleurs par leur direction d'arrivee, le temps, le spectre ou la puissance 
comme cela est decrit dans Touvrage de R.A. Monzingo. T.W. Miller, intitule "Introduction to adaptive 
arrays", edite par John Wiley and Sons, New-York, 1 980 et la these de P. Chevalier, intitulee "Antenne 
adaptatrve : d'une structure lineaire a une structure non lineaire de Volterra", These de doctorat, Uni- 
versHe de Paris Slid, Juin 1 991 puis de correler le signal de reference avec la sortie de I'antenne. Les 
35 techniques de f iltrage spatial discriminant le signal et les brouilleurs par le temps, le spectre ou la puis- 

sance sont generalement sous-optimales, alors que celles exploitant la connaissance de la direction 
d'arrivee du signal necessitent la mise en oeuvre d'une technique de goniomfetrie, coflteuse et delicate 
dans la gamme des frequences HF notamment du fait des multi-trajets. 
Pour ces raisons, les techniques de la seconde famille. visant a effectuer la rejection de brouilleurs conjoin- 
40 tement avec la prise de synchronisation, a partir de la connaissance du signal de reference, appara.ssent plus 
prometteuses Une telle technique a par exemple ete proposee par M. Compton dans son article "An adaptive 
array in a spread spectrum communication system" paru dans Proc. IEEE. vol. 66, pp. 289-298, March 1978, 
pour la synchronisation de signaux PN, a partir d'une architecture analogique integrant des "Delay Lock Loop". 
Une seconde technique fut proposee par L.E. Brennan et I.S. Reed dans leur article intitule "An adaptive 
45 array signal processing algorithm for communications", paru dans IEEE Trans-Aero-Syst. vol. AES-18. n°1. 
pp 1 24-1 30. June 1 982 pour les transmissions numeriques et completee par D.M. Duglos et R.A. Scholtz dans 
I'article intitule "Acquisition of spread spectrum signals by an adaptive array", paru dans IEEE Trans-Acous- 
Speech-Signal-Proc, vol. ASSP-37. n°8. pp. 1253-1270. August 1989. Ces techniques apparaissent toutefois 
beaucoup plus couteuses que les.techniques mono-capteur et sont malheureusement inapplicables telles quel- 
50 les en presence des derives Doppler du signal inherentes aux transmissions ionospheriques. 
Le butde I'invention est de pallier les inconvenients precites. 

A cet effet, I'invention a pour objet un procede permettant a un modem de se synchroniser sur un trans- 
metteur de signaux numeriques par voie hertzienne en presence de brouilleurs consistent a: 

- estimer la matrice de correlation Rxx de signaux recus sur un ensemble de N capteurs de reception par 
55 blocs composes d'un nombre d'echantillons d'une sequence de synchronisation transmise par le trans- 
metteur, A 

- a calculer la matrice de correlation inverse R«, 
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- a calculer des vecteurs d'intercorrelation f Xd entre les signaux X(k) recus sur I'ensemble des N capteurs 
et un signal de replique d(k) connu, 

- a calculer un critere de synchronisation ef fectuant le produit scalaire entre le vecteur d'intercorrelation 
rxd et le vecteur obtenu par produit de I'inverse de la mat rice de correlation et du vecteur d'inter- 

5 correlation r Xd - . 

- a comparer ia valeur du critere obtenu a une valeur de seuil t| determinee pour placer la synchronisation 
sur I'echantillon pour lequel la valeur des criteres depasse la valeur de seuil j\ t caracterise en ce qu'il 
consiste pour estimer le d6calage Doppler inherent aux communications entre emetteurs et recepteurs, 
a decouper la sequence de synchronisation en un nombre determine de M tranches de ko symboles, a 

10 effectuer sur ces tranches des calculs de vecteurs_ d'intercorrelation partielles r Xd i (i ^ 1, M) de telle 

sorte que sur chacune des tranches de ko symboles, le decalage Doppler fasse tourner la phase d'une 
valeur inferieure a 90°, a effectuer des calculs de transformee de Fourier rapide (FFT) sur les vecteurs 
?xdi obtenus, et a calculer une valeur du critere de synchronisation a partir des vecteurs r Xd {kAf), obtenus 
en resuitat des calculs de trarisform6e de Fourier rapide, et de la matrice de correlation Rxx* 

15 L'invention a egalement pour objet un dispositif pour la mise en oeuvre du procede. 

Le procede et le dispositif selon l'invention ont pour avantages qu'Ms permettent a la fois d'assurer la syn- 

- chronisation d'un modem dans des conditions de brpu illage, severe. et en pr6senpe de multi-trajets, d'estimer 
*-*"* la derive en frequence entre emetteuret rebepteuretde detecter les instants d'arrivee detous les multi-trajets 

— /assoctes au signal utile, pour vu que ceux-ci.soient de puissanpe suff isamment forte par rapport a la puissance 
20 de bruit de fond. : . 

D'autres caracteristiques et avantages de ('invention apparaTtront dans la description qui suit faite au re- 
gard des figures des dessins annexes qui repr^sentent : 

- La figure 1 un schema fonctionnel d'une antenne adaptative selon l'invention. 

- La f igure 2 un exemple de mise en oeuvre d'un filtre adapte spatial a discrimination par la modulation 
25 pour la realisation d'une antenne adaptative selon rinvention. 

- La figure 3 un schema d'antenne adaptative pour lajnise en oeuvre du procede de synchronisation spa- 
tiale selon rinvention. 

- La figure 4 un graphe pour illustrer la methoded'estirnation de la matrice de correlation de signaux pro- 
venant de N capteurs. - 

30 . |_ a figure 5 un graphe pour illustrer le d6coupage des symboles de la sequence de synchronisation par 

tranche pour le calcul de vecteur d'intercorrelation en presence de d6calage Doppler. 

- La figure 6 un graphe montrant Involution de critere de synchronisation mis en oeuvre par l'invention 
en presence d'un decalage Doppler. 

- La figure 7 un organigramme retracant Tensemble des etapes necessaires a la mise en oeuvre du pro- 
35 C ede de synchronisation selon l'invention. 

Le precede de synchronisation mis en oeuvre par l'invention fait appel aux techniques connues de syn- 
chronisation mono-capteur pour les appliquer a la realisation d'antennes adaptatives sur replique. 

Selon la technique de synchronisation mono-capteur, le critere de synchronisation C(n) est obtenu en cal- 
culant la correlation y(n) entre K symboles connus dans une sequence de synchronisation et les K echantillons 
40 qui suivent I'arrivee d'un echantillon n sur le capteur, en normalisant les r6sultats par la puissance p(n) du signal 
recu. Ces calculs sont effectues en appliquant les relations : 

45 y(n) = lf>(k + n)d-(k) (1) . 

K k=i 



50 



55 



1 K 

p(n) = ^Xx(k+-n)x*(k + n) (2) 
K k=i 

ou: 

- x(k) est le signal issu du capteur utilise, 

- d(k) est le signal de replique connu (K symboles places en debut de trame), 
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■ ' ■ *- * represente Toperation de conjugaison sur lesnombres complexes.; 

Ces calculs sont reiteres pour chaque echantillon n, et lorsque le entire C(n) depasse un certain seuil i\, 
dependant de la probability de fausse alarme, la synchronisation est decidee. L'instant de synchronisation est 
determine parte rang d'echantillon n 0 pour lequel C(n 0 )> i\. La normalisation par la puissance p(n) permet de 
5 choisir un seuil independent du niveau de signal recu. 

Dans ce procede la synchronisation est obtenue lorsque le rapport signal a bruit est suffisant. Ses per- 
formances dependent en partie du nombre d'echantillons pris en compte pour la correlation. Cependant en 
presence de brouilleurs forts, la synchronisation ne peut etre effectuee. 

L'idee mise en oeuvre dans I'invention consiste alors a effectuer la synchronisation, en sortie d'un traite- 
10 ment d'antenne adaptative, permettant la rejection des brouilleurs. 

Le schema fonctionnel d'urie antenne adaptative est represente a la figure 1 , celle-ci comporte un reseau 
1 de N capteurs 1 ^ ... 1 N , un reseau de filtres adaptes 2 (constitues par des circuits de ponderations amplitude- 
phase, des filtres a reponse impulsionnelle f inie RIF ou inf inie Rll) couple a un processeur de commande adap- 
tatif 3.' Le processeur 3 calcule les coefficients W 1 a W N du reseau de filtres 2 pour I'adapter au signal utile 
15 et aux brouilleurs, en I'absence de toute connaissance a priori sur les brouilleurs et a partir d'une information 
minimale sur le signal utile. Les coefficients Wi a W N permettent de former un trou dans le diagramme de 
rayon nerhent du reseau d'antenne dans la direction des brouilleurs, tout en maintenant un certain gain dans 
la directibhdu signal utile. -v-^ ' 

Differents types d'iriformatiori sur ^ le signal utile sont exploites par rantenne adaptative. Ceux-ci sont d'or- 
20 dre spatial lorsque la direction d'arrivee du signal utile est connue, d'ordre temporelle lorsque les instants de 
presence et d'absence de signal utile sont connus par exemple en evasion de frequence, relatifs a la forme 
d'onde lorsqu'une information sur la modulation du signal utile est disponible pour disposer d'une replique du 
signal utile en intercalant par exemple des sequences connues dans la trame utile, d'ordre spectral lorsque 
le signal utile occupe une certaine bande de frequence et que les brouilleurs occupentune bande plus impor- 
25 tante, ou statistiques pour separer les differentes sources arrivant sur le reseau en exploitant par exemple les 
statistiques d'ordre 4 des signaux en supposant que les differents signaux sont statistiquement independants. 

Dans le cas de Tantenne adaptative representee a la figure 2, P information sur le signal utile est constituee 
par les sequences de symboles connus intercalees dans la trame entre les sequences de symboles d'infor- 
mation. Les filtres places en sortie de chacun des capteurs 1 , a 1 N sont formes respectivement par des circuits 
30 de ponderation amplitude-phase \N, a W N des signaux Xl (t) a x n (t) fournis par les capteurs. Les signaux obtenus 
en sortie des circuits de ponderation sont appliques a un circuit de sommation 4. 

Les caracteristiques du reseau de filtres 2 sont calculees en minimisant une Erreur Quadratique Moyenne 
(EQM) estimee £ entre un signal de replique d(k) et le signal y(k) sortant de I'antenne de la facon schematisee 
^ la figure 2 ou les elements homologues a ceux de la figure 1 sont reperes avec les memes references, § est 
35 definie par la relation 



K K o 

1 = Z|y(k) - <k)| 2 = Z W*X(k) - d(k)f (4) 
k=1 k=1 



ou: 



y(k) est la sortie de I'antenne adaptative et est definie par 

y(k) = WX(k) (5) 

45 -W est le vecteur poids def inissant I'antenne adaptative, W 1 etant le vecteur transpose conjugue corres- 

pondant. (le signe t en exposant def init Toperateur correspond ant.) 

- XCkj'est le vecteur forme par les signaux recus sur les capteurs. 

- d(k) est le signal de replique connu. ^ 

La difference y(k) - d(k) est obtenue sur la figure 2 a la sortie d'un circuit soustracteur 5. 
so Le vecteur poids W est calcule par la formule connue de Wiener : 

W = R-L?x d (6) 



avec : 



RxX = ^ZX(k).X t (k) (7) 
K k=1 
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matrice de correlation des signaux provenant des capteurs 1 , a 1 N . (X*(k) etant le vecteur transpose conjugue 
du vecteur X(k)). * 



45 



? xd=^i;x(k) d(k) 

K k=l 



vecteur d'intercorrelation entre X(t) et d(t) (8). 

10 L'antenne adaptative ainsi utilisee conduit a une maximisation du rapport signal/(bruit + interference) 

connu sous la designation SNIR abreviation anglo-saxonne.de ^Signal to Noise + Interference Ratio" obtenu 
en sortie de l'antenne, avec une possibility de rejection pouvant aller jusqu'a N-1 brouilleurs, et est appelee 
Filtre Adapte Spatial (FAS). 

La recherche de la synchronisation est effectuee selon Tinvention en sortie de l'antenne adaptative en 

15 discriminant le signal selon sa modulation, et en utilisant le critere de synchronisation mono-capteur. Le sche- 
ma de principe de la synchronisation spatiale est represente a la figure 3 ou les elements homologues a ceux 
des figures 1 et 2sont : represents avec les memes references^ Le critere de synchranisation.cieyient en uti- 
lisant les notations precedentes : : 



55 



C(n) = l|v(k + n)d*(k) (9) 



K k=l 



25 y(k+n) represente la sortie de I'antenne adaptative conduisant a une minimisation de I'erreur quadratique 

moyenne L(n) estimee sur K echantillons entre la replique d(k) et le signal arrivant pris a partir de I'echantillon 



moyenne £(n) 
n, suivant la relation 



E(n)= ^|y(k n) - d(k)| 2 avec y(k + n) =.W(n) t .X(k + n) (10) 
k=1 



avec : A 
35 W(n) = Rxx- 1 (n).r X d(n) (11) 

Les differentes correlations sont calculees sur les K echantillons qui suivent Tarrivee d'un echantillon n 
suivant les relations : 



Rxx(n) = ^ ZX(k + n)X l (k + n) (12) 
K k=l 



Wn) = ^£X(k + n)d*(k) (13) 
K k=i 



Lorsque I'echantillon n correspond a la position de synchronisation, les poids W(n) def inissent le filtre 
so adapte, et y(k+n) contient alors le signal utile debarrasse de la contribution des brouilleurs. Dans le cas contrai- 
re, aucun signal n'est correle avec ia replique d(k) et y(k+n) ne contient plus alors que la contribution du bruit 
de fond. 

L'expression du critere de synchronisation spatiale se simplifie en ecrivant : 



K ( K "\ 

C(n) = -1 XW(n) t X(k+n)d*(k) = W(n) t ± £x(k + n) d* (k) 
K k=1 l K k=1 ) 



(14) 
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soit A 

C(n) = ^ d (n)Rxx(n)- 1 f Xd (n) (15) 

De meme que pourle casmono-capteur, ce calcul est reitere pourchaque echantillon n, etiorsque le critere 
C(n) depasse un certain seuil n, dependant de la probabilite de fausse alanme, la synchronisation est decidee 
(la normalisation par la puissance est ici obtenue par rintermediaire de R^ 1 ). 

L'instant de synchronisation est alors determine par le rang de Techantillon pour lequet C(n 0 ) est superieur 
a une valeur de seuil determinee tj. 

Les calculs correspondant sont effectues sur la figure 3 par le bloc 7 de calcul du critere C(n) et par le 
detecteur de seuil 8, a partir de la sequence de synchronisation d(t) et du signal y(t) correles par le correlateur 
6. 

Cela necessite pour chaque echantillon, d'une part de calculer N correlations entre le signal de replique 
d(k) et les signaux x(k) repus paries N capteurs (sort N fois les calculs effectues pour la synchronisation mono- 
capteur). et d'autre part, d'estimer puis d'inverserla matrice d'intercorrelation des signaux capteurs. 

Au cours de ces calculs la synchronisation et le calcul du f iltre adapte 2 sont effectues conjointement. 
15 Lorsque la valeur du critere depasse le seuil t|, calcule, la synchronisation est decidee et le vecteur poids de- 
finissant le f iltre 2 s'en deduit facilement par la fprmule connue de Wiener. Le vecteur poids W peut ensuite 
etre utilise pour filtrer les 6chantillons contenant I'information. . ; : ■ . ; 

,En supposant que la replique d(t) est de puissance ijes correlations tendent vers les valeurs suivante^ 

iorsque le nombre d'echantilibns K devient tres grand : : 71 r * 
20 a ,r~ c 

ou: 

- R est la matrice de correlation des interferences et du bruit de fond, 
25 - S est la matrice du vecteur directeur du signal utile, 

- 7t s est la puissance du signal utile regu sur un capteur. 

Le critere de synchronisation qui converge alors vers une valeur C 0 s'ecrit encore: 

S*R~ 1 S = SNIRopt {16 , 
0 1 + S^R- 1 S 1 + SNIRopt V ' 
30 ou SNIRopt est le SNIR asymptotique obtenu en sortie du f iltre spatial adapte 2. 

Lorsque le SNIR est fort, le critere tend vers 1 , lorsquMI est faible, il tend vers 0. L'importance du choix du 
seuil r|, fixe la probabilite de detection de la synchronisation. Le seuil depend naturellement de la probabilite 
de fausse alanme souhaitee. A 

Le calcul du critere de synchronisation C, de ? Xd et de Rxx est effectue pour chaque position de la syn- 
35 chronisation et notamment la matrice Rxx est estimee et inversee egalement pour chaque position de la syn- 
chronisation, la puissance de calcul demandee est elevee. Mais il est possible de reduire la puissance de calcul 
en estimant la matrice de correlation Rxx moins souvent sur un intervalle de temps plus important (K' 6chan- 
tillons au lieu de K avec K>K de la fa$on representee a )a figure 4) et en utilisant cette estimation de Rxx sur 
toutes les positions de synchronisation incluses dans cet intervalle de temps. 
40 Dans ces conditions la matrice Rxx s'ecrit : 



45 



2B(t k )B l (tk)+ Ss^)^) 



(17) 



k=1 



k=K n 



so ou : 

- B(t k ) est la contribution du bruit et des brouilleurs sur X(t k ), 

- S(t k ) est la contribution du signal utile sur X(y, 

- Ko est I'instant d'apparition du signal utile sur le palier, 

- e est un terme prenant en compte Tintercorrelation entre le signal B(t) et le signal utile S(t), qui tend vers 
55 0 lorsque K devient grand. 

Pour la position de synchronisation testee : 
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fXd = K 



K 0 4K.. 

Z X(tk)cf(k) 

k=Ko 



(18) 



10 



15 



20 
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Dans les conditions asymptotiques (avec K»1), lorsque te signal utile est synchronise avec la repiique 
d(t), Rxx et ^xd se simplif ient et s'ecrivent : 

R XX ^R + K> -^ + 1 tt s SS* (19) , / 

fxd^V^S (20) 

Le critere obtenu en estimant la matrice de correlation Rxx sur K* echantillons possede comme valeur 
asymptotique, apres un calcul utilisant le iemme d'invefsion matricielle: 



C'n = 



SNIRopt 



1 ^ K'-Ko 1 SN|RQpt 



(21) 



Comme * 1 < 1y alors C'b > C 0 K • ' r :^7 : " " -V ' 1 ' : , > :. ff . ; , - 

La valeur du critere apparaft dans ces conditions augmentee d'un faeteur .dependant de Tiristant d'appa- 
rition du signal utile sur la frequence consideree, ce qui entrame une amelioration de la prise de synchronisa- 
tion, au prix d'une puissance de calcul reduite. 

D'autre part, le nombre K' ne doit pas etre choisi trap grand, sa valeur est determinee pour que sur 
K' Echantillons, I'environnement de brouillage puisse etre considere comme stationnaire. 

Le seuil a prendre en consideration depend de la probabilite de fausse alarme (Pfa). La (Pfa) est la pro- 
babilite qu'un bruit soit detecte comme une synchronisation. Le calcul de la valeur du seuil r\ a lieu en calculant 
Pfa=P(C>ri) lorsque I'antenne recoit uniquement du bruit. 

En presence de bruit le critere de synchronisation s'ecrit: 



30 



K 



^, K 



k=1 



^SBkBkM.UZBk^j (22) 



En supposant que B k est un bruit gaussien complexe de puissance o 2 , independant entre les capteurs, 
^ 5 lorsque K est suff isamment grand, le terme du milieu de la relation (22) devient : 



40 



jlBkBu'k^M (23) 



(Id etant la matrice identite on unitaire) 

le critere de synchronisation peut alors etre approxime de la facon suivante : 



45 



K 2 a 2 



( K 



( K 



(24) 
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Soit : 
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N 



c = 



K a i=1 



K 

Z B k < 

k=1 



|2^ 



s ^Zfo 2 + * 2 ) ( 25 ) 



K a i=1 



7 



avec: 



P 
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K 



Z = R, + j Ij (26) 

k=1 



Bjc represente le bruit obtenu sur le capteur i, pour I'echantillon k. 

Ri + jl| peut etre approxim6 par un bruit gaussien complexe de variance K a 2 . Dans ces conditions Rj et I, 
1° sont chacun des bruits gaussiens reels de variance 



15 



En notant T, = Rj cu l jf C devient : ....... 



20 



2N 



K 



2 Jl 



Z(li 2 ) (27). 



i-1 
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K <T 2 

ou T| est une variable aleatoire qui suit une loi normale centree de variance — — . La variable aleatoire Z f = 



<2 



suit done une loi normale centree de variance 1. 



Wo 2 

C devient : 



- 2N 

Cs 5cS2f (28) 



1=1 



35 



Comme la suite des Zj peut etre consideree comme une suite de 2N variables al6atoires independantes 
gaussiennes de moyennes nulles et de variances 1, la variable aleatoire 



40 



2N 

i=1 



45 



suit une loi du x 2 & 2N degr6s de liberty. 
La probability d'avoir C ^ n s'Scrit alors : 



p(c > n) = p 



50 



(2N - 
£zf>2K n 
0=1 



= Pfa (29) 



Le tableau suivant donne differentes valeurs de seuils t| obtenues par recherche dans une table donnant 
P(Z^A) pour une loi du x 2 a 2N degrds de liberte et pour Pfa = 10" 5 . 
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N = 5 
N = 4 
N = 3 
N = 2 
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A = 42 
A = 37 
A = 33 
A = 28 



t! = 0,26 
t, = 0,23 
ti = 0,21 
^ = 0,18 



Comme entre remission et la reception, le signal transmis peut etre affecte d"un certain decalage en fre- 
quence, ou decalage Doppler, du essentiellement aux derives en frequences des oscillateurs locaux des 
10 synthetiseurs a remission et a la reception, ainsi qu'aux conditions de propagation ionospherique, il est ne- 
cessaire de compenser ce decalage. Le risque est s'il n'est pas compense. d'empecher toute prise de syn- 
chronisation. 

En appelant Af 0 , le decalage Doppler entre le signal regu parl'antenne etle signal attendu eten considerant 
la position de synchronisation, le vecteur signal recu peut s'ecrire avec les notations des paragraphes sui- 
ts vants : 

X(t) = d(t) e 2 )"^ S + B(t) 
- " -La -matrice de. oorrelatipn estimee peut alors.*'$crire •: . O /.v'-- *■.■<<& 



20 



RXX - £ Z <t k )e 2 j" Af o t K s i B(tK)|^ d(t k )e 2 -j- Af o t K.S + B^)]* (30) 
k = 1 
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K 

^XX = « l + if Z^s|d(tk)| 2 S S 1 + B(t k ) B^) 1 "] (31) 
K k=1 L J 



ou e est une matrice prenant en compte les interconrelations entre la contribution du signal utile et le bruit b(t). 
30 Lorsque K»1, dans les conditions asynrvptotiques, Rxx peut s'approximer par : 

Rxx = «s S S + + <r 2 Id (32) 
La relation (32) montre que le decalage Doppler ne modif ie pas le calcul de Rxx- En effet, le decalage 
s'elimine nature II erne nt lorsque le signal utile est multiplie par sa quantite conjuguee. 
Le vecteur d'intercorrelation peut s'ecrire : 

35 



p Xd = ^ Z[V^ d^**-** S + B(t k )] cf (tk) (33) 



? Xd = e' + 1 ij^Kt^f eW^-tk S ] (34) 



ou e' est un terme prenant en compte I'intercorreiation entre le bruit B(t) et la replique. Lorsque K»1, r Xd peut 
s'approximer par : 



50 



Af 0 



K k=l 



S (35) 



55 La relation (34) montre que ie decalage Doppler modifie I'estimation du vecteur d'intercorrelation r Xd et 

que du fait que la phase du terme en exponentiel complexe peut tourner de plusieurs fois 360° sur I'intervalle 
de temps utilise pour le calcul de la correlation. ('integration ne peut plus avoir lieu de fagon coherente ce qui 
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empeche le critere de synchronisation multicapteurs de converger. Le processus suivant permet d'estimer le 
decaiage Doppler et de le compenser pour rendre possible la prise de synchronisation. 

Selon ce processus, les K symboles de la sequence de synchronisation sontdecoupes de la facon repre- 
sentee a la figure 5 en M tranches de Ko symboles (M Kq = K), de telle sorte que sur chacune des tranches 
de Ko symboles, le decaiage Doppler fasse tourner la phase du signal regu par Tantenne d'une valeur inferieure 
a 90°, af in de permettre une integration coherente de rxd sur chacune des tranches. 

Si r X di designe la valeur de i'intercorrelation estimee sur la tranche i, celie-ci s'ecrit : 



10 



r Xdi = 



* Z J^kjfe F < 
^ k=iK 0 +1 



S i = 0....,M-1 (36) 



15 



En mettant en facteur dans la relation (36) un terme dependant uniquement de Ko, la longueur de la tranche 
elementaire considered, et du decaiage Doppler, la relation (36) devient : * ^ 



20 



r Xdi = 



^Z^K(^iKo)T)| 2 e F * 



2j7t.^.iK 0 



e Fe .S i = 0,„.,M-1 (37) 



25 



30 



Ce qui peut encore s'ecrire 



fxdi ^ X|-e3*-ir-«<o .S i = 0 M - 1 (38) 



La relation (38) montre que les vecteurs r X di sont des vecteurs, echantillonnes a Leur transformee de 
Fourier discrete est alors : 



35 



M-1 



f Xd (Af) b X r Xdi e * Fe/K o\ (39) 
i=0 



40 



M-1 2j7c- f °~ Af: 



0Cd(Af)s £ Xje Fe 
i=0 



iKo 



(40) 



La relation (40) montre que lorsque Af = Af 0 , le decaiage Doppler est compense, ce qui rend maximum la 

valeur du critere de synchronisation r Xd (Af) t Rxx" 1 rxd(Af)- 

Af in d'estimer Af, c'est-a-dire af in d'estimer le decaiage Doppler qui conduit a la valeur maximale du cri- 
^ tere, et af in de reduire le cout de calcul correspond ant, une transformee de Fourier rapide est effectuee sur 
les vecteurs rxdi (le nombre de tranches elementaires M etant choisi pour former une puissance de 2) : ce calcul 
donne : 



50 



M— 1 — «.jjt. .in, ivi-i 

r^kAf^ £ r X die = £ f Xdi e k = 0 M-1 (41) 

i=0 



-2)71, 



-.ik M-1 

z 

i=0 



0 - ik 



oil Afi est la resolution Doppler telle que : 

<«> 

Si la resolution est insuff isante, une amelioration d'un facteur 2 peut etre obtenue en calculant une trans- 
formee de Fourier rapide, sur les vecteurs r Xd i "decalees" en frequence de la facon suivante : 



10 



-ER0 665 665 A1 



r Xd kAf, + 



M-1 

Z 

i=0 



-J7l.— 



e M 



k = 0,...,M-1 - (43) 



L'operation consiste a multiplier chacun des vecteurs d'intercorrelation r Xd s par le terme e- 0*^) et a effec- 
tuer ensuite un calcul de transformee rapide sur les vecteurs ainsi obtenus. 

Le critere de synchronisation est ensuite estime sur chacune des cases de calcul des transformees de 
w Fourier precedentes par la relation : 



15 



Bien eyidemment, sj la resolution est toujours insuffisante, il est possible de rameliorer en utilisant le 
meme processus pour calGuleir des vecteurs 
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f Xd (kAf,-3) 
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45 
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ou n est un entier quelconque. 

Un test sur la case de la transformee correspondant au maximum du critere de synchronisation est effec- 

tue. Si le seuil est depasse pour une valeur ko donnee du decalage Doppler, il est alors possible d'aff iner 

le calcul du decalage Doppler (indispensable pour la demodulation) en calculant la valeur du critere pour dif- 

ferentes valeurs de decalage Doppler entourant ko En effet, a Tissue de cette etape, la resolution sur I'es- 

timation du decalage Doppler est Pour Tameliorer, il faut proceder par dichotomie : a chaque iteration et 

de la facon representee a la figure 6, la valeur du critere pour la valeur de decalage Doppler situee au milieu 
des deux valeurs de decalage Doppler donnant les valeurs maxima du critere lors de Iteration precedente 
est calculee. La recherche par dichotomie est arretee des que la resolution souhaitee sur I'estimation du de- 
calage Doppler est obtenue. Apres n iterations de I'algorithme de recherche par dichotomie, la resolution, qui 



etait initialement de est : 



AM 2 = Af-| 
2" 2 n + 1 



(45) 



Uenchamement des etapes de calcul precedente pour la mise en oeuvre du precede selon Tinvention a 
I'aide de moyens de traitement appropries formes par exemple par des processeurs de traitement du signal 
ou tout dispositif equivalent connu est figure sur Torganigramme de \a figure 7. 

II consiste suivant Tetape 9 a estimer la matrice de correlation Rxx suivant la relation (7). Cette matrice 
est estimee sur un bloc contenant un nombre d'echantillons K' superieur au nombre d'echantillons de la se- 
quence de, synchronisation. 

^estimation de Rxx est utilisee ensuite pour toutes les positions de synchronisation incluses dans le bloc. 
Puis la matrice de correlation inverse R^ 1 est calculee a I'etape 1 0. r Xdi qui est def ini par la relation (8) est cal- 
cule a I'etape 11. Pour chaque echantillon correspondant a une position de synchronisation P dans la suite 
des echantillons P = 1 a P = K'-K, f X di est calcuie pour i = 1 a M sur M blocs de Kq symboles comme defini a 
la figure 5 suivant la relation 
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1 



Kn 



rxdi = tt Z x ( p +K + iK o) d (K + iK 0 ) (46) 
K o K=1 



11 



niietrw^irv. -cn neeceec* i I - 
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Des calculs de transformee de Fourier Rapide sont effectues £ retape' suivante 12 sur un nombre deter- 
mine M (M = 8 par exemple) de vecteurs f Xd i suivant la relation (41). 

Aux etapes 1 3 et 14 un calcul de vecteuir R Xdi decal6 et un calcul de transformee de Fourier rapide suivant 
la relation (43) sont executes pour ameliorer la resolution par un facteur de 2. 

Le critere de synchronisation est calculi d I'etape 1 5 en application de la relation (44). 

La comparaison du entire og seuil t| est effectuee a I'etape 16. Si & I'etape 16 la valeur du critere 

' <m 

est superieure ou egale au seuil, la position de synchronisation est definie d I'etape 17 sur la position P de 
I'echantillon correspondant ^estimation du decalage Doppler peut ensuite etre aff ine & I'etape 1 8 en calculant 
a nouveau la valeur du critere de synchronisation pour des decalages Doppler entourant le decalage 

precedent. Cette etape peut etre suivie par une etape 19 d'estimation de la position de synchronisation pour 
les autres trajets. 

La recherche de synchronisation prec6demment decrite est effectuee pour chacun des echantillons de 
signal numerise. En pratique, pour ameliorer la precision de la position de la synchronisation il est preferable 
d'effectuer un surechantillonnage du signal d'un facteur determine L, L etant egal & 4 par exemple. Les trai- 
tements de I'organigramme de ia figure 7 et decrits ci-dessus sont encore utilisables sur le signal ainsi sure- 
chantillonne, mais il est possible d'utiliser une variante dans le precede, af in de reduire la puissance de calcul 
necessaire au traitement. Dans ces conditions, apres avoir calcuie pour chaque echantillon le vecteur r Xd i pour 
i = 1, ... M sur M blocs de Ko symboles suivant la relation (46) un critere de synchronisation Ci est calcuie sur 
chacun des M blocs par la relation 

Ci = r\ dr R^.f Xdl (47) 
Puis un critere spmme de Tensemble des Cj est calcuie suivant la relation : 

Cl(P) = ZCi (48) 

. _ i=i 

Ce critere correspond & une integration non coherente des Ci (contrairement & ce qui est fait pour le calcul 
par FFT ou on integre de fagon coherente), ce qui se traduit par une remontee du critere sur les echantillons 
de bruit. 

Le calcul du critere a lieu par blocs de L echantillons consecutifs. Sur chacun de ces blocs, I'echantillon 
qui donne la valeur maximum de C^P) est selectionne pour effectuer le calcul par transformee de Fourier ra- 
pide necessaire a Testimation du decalage Doppler comme deja decrit aux etapes 11a 19 de la figure 7. 



Revendications 

1. . Precede permettant a un modem de se synchroniser sur un transmetteur de signaux numeriques parvoie 
hertzienne en presence de brouilleurs consistant & : 

- estimer (9) la matrice de correlation R^ de signaux regus sur un ensemble de N capteurs de recep- 
tion (1i ... 1 N ) par blocs composes d'un nombre d'echantillons d'une sequence de synchronisation 
transmise par le transmetteur, A 

- a calculer (10) la matrice de correlation inverse Rxx\ 

- a calculer (11) des vecteurs d'intercorrelation f Xd entre les signaux X(k) re$us sur ('ensemble des N 
capteurs et un signal de replique d(k) connu, 

- & calculer un critere de synchronisation (7) ef fectuant le produit scalaire entre le vecteur d'intercorre- 
lation r xd et le vecteur obtenu par produit de ('inverse de la matrice de correlation R^ et du vecteur 
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d'interconr6latipn.rxd 

- a comparer (8) la valeur du critere obtenu a une valeur de seuil tj determined pour placer la synchro- 
nisation sur I'echantiilon pour lequel la valeur du critere depasse la valeur de seuil r\, caracterise en 
ce qu'il consiste pour estimer le decalage Doppler inherent aux communications entre emetteurs et 
recepteurs, a decouper la sequence de synchronisation en un nombre determine de M tranches de 
ko symboles, a effectuer (11, 13) sur ces tranches des calculs de vecteurs d'intercorr6lation partielles 

f Xdi (j = 1 M) de telle sorte que sur chacune des tranches de ko symboles, le decalage Doppler 

fasse tourner la phase d'une valeur inferieure a 90°, a effectuer (1 2, 1 4) des calculs de transformee 
de Fourier rapide (FFT) sur les vecteurs f X di obtenus, et a calculer une valeur du critere de synchro- 
nisation a partir des vecteurs r Xd (kAf) A obtenus en r6sultat des calculs de transformee de Fourier ra- 
pide, et de la matrice de correlation Rxx- 

A 

Precede selon la revendication 1 caracterise en ce que Testimatidn de la matrice de correlation Rxx est 
obtenu par la relation 

. ■ Rxx =4S x(k) x w 

et le vecteur d'intercorrelation est estime suivant la relation 



1 

r Xd = ^Z X (k + p) d ( k )/ 
^ k-1 

ou p est la position du premier 6chantillon a partir duquel la correlation est effectuee entre le vecteur X 
et la replique d. 

Procede selon la revendication 2 characterise en ce qu'il consiste a choisir la case de calcul de la trans- 
formee rapide qui donne une valeur maximum du critere de synchronisation et a comparer la valeur maxi- 
mum obtenue a une valeur de seuil determinee pour decider de la synchronisation. 

Proc6d6 selon Tune quelconque des revendications 1 a 3 caracterise en ce qu'il consiste pour obtenir une 
meilleure resolution sur le decalage Doppler et une valeur plus importante du maximum du critere, a cal- 
culer (13) des vecteurs r xdi decaies en frequence et a effectuer un calcul de FFT sur les vecteurs' r xdi de- 
cales de facon a obtenir la valeur du critere avec une meilleure resolution sur les frequences Doppler. 

Proc6d6 selon Tune quelconque des revendications 1 a 4 caracteris6 en ce qu'i I consiste a sur6chantil- 
lonner les signaux regus sur les N capteurs par rapport au rythme symbole, pour ameiiorer la resolution 
sur la position de synchronisation. 

Precede selon Tune quelconque des revendications 1 a 5 caracterisee en ce qu'il consiste a effectuer le 
calcul de la matrice de correlation Rxx sur un nombre d'echantillons K* superieur au nombre K d'echan- 
tillons de la sequence de synchronisation af in de diminuer la puissance de calcul. 

Procede selon Tune quelconque des revendications 1 a 6 caracterise en ce qu'il consiste pour diminuer 
la puissance de calcul et apres avoir surechantillonne d'un facteur L par rapport au rythme symbole, a 
effectuer une etapede pretraitement consistant a calculer pour chaque echantillon p(p = 1 aL)composant 
un symbole les valeurs C ( du critere de synchronisation sur chacune des M tranches de ko symboles et 
a calculer la valeur du critere en effectuant la somme des criteres Q obtenus pour retenir dans la suite 
du traitement rechantillon donnant la valeur maximale du critere C = ZC t . 

Dispositif pour la mise en oeuvre du procede selon Tune quelconque des revendications 1 a 7 du type 
comportant, un r6seau de N capteurs (1 ,...1 N ) de signaux de reception couple a un reseau de f iltres adap- 
tes (2) comportant des circuits de pond6rations amplitude phase commandes par un processeur de 
commande (3), le processeur de commande (3) etant programme pour 
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- estimer la mat rice de correlation inverse R^J des signaux recu X (K) et calculer les vecteurs d'inter- . 
correlation r Xd entre les signaux X(k> recus et le signal de replique d (K) caracterise eri ce qu'il comprend 
des moyens de traitement (5.„19) pour calculer un critere de synchronisation (7) en effectuant un 
produit scalaire entre le vecteur d'intercorrelation r xd et le vecteur obtenu par le produit de Tinverse 

5 de matrice de correlation R^J et du vecteur d'intercorrelation r^, 

- comparer (8) la valeur du critere obtenu a la valeur de seuil r\ pour placer la synchronisation sur 
Techantillon pour lequel la valeur du critere depasse la valeur du seuil r\, 

- et estimer le decalage Doppler inherent aux communications entre emetteurs et recepteurs, en de- 
coupant la sequence de synchronisation en un nombre determine de M tranches de ko symboles, 

10 en effectuant (11, 13) sur ces tranches des calculs de vecteurs d'intercorrelation partielles f xdj (i = 

1, ...,M) de telle sorte que sur chacune des tranches de Kq symboles le decalage Doppler fasse tour- 
ner a la phase d'une valeur inferieure a 90°, en effectuant (12, 14) des calculs de transformee de 
Fourier rapide (FFT) sur les vecteurs f Xd , pbtenus, et en calculant une valeur du critere de synchro- 

15 nisation a partir des vecteurs r X d(KAf[, obtenu s en resultat des calculs de transformee de Fourier ra- 

pide et de la matrice de correlation Rxx- 

> * ■ i ' ■' . 
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Calcul du critere 



^7 



CIn)4^ y(k + n)d*(k] 
C[n) = r* d (n)R; 1 x (n)f Xd (n) 



Synchronisation- si •. 

C{n) >rj 
lr\ dependant de la Probabilite 
de fausse alarme desiree) 



FIG.3 
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Estimation de R xx sur K' echantillons 



K echantillons* 



K echanti lions 



1«J« position tester avec R xx 
estimee sur K* echanti lions 



Darniere position^ testee avec R xx 
estimee sur K 1 echantillons 



K 0 



xdl 



FIG. 4 
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FIG.5 
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valeur de decalage Obppler estimee 

FIG. 6 
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Estimation de. R xx 
sur K' echantillons 




non 



Fin 



Calcul de- R 
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P=1 

3: 



Calcul' .de' f XC ji^ 
i= 1, ....... M- 2* 
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Calcul de FFT sur 
tes vecteurs f xt jj 
' ^dtk^oj 

~l 
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Calcul des vecteurs 
r xdj decales 

1 



Calcul de FFT sur ^ 
les vecteurs r xd j de- 
cales ^r xd (k^o+Afo) 



i 
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Calcul du critere C (k^- f o) 
k =1, , 2M 
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2 * 
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P>K'-K 




P= P+1 


^ non 
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